High- resolution detection of Brownian motion for quantitative Optical Tweezers experiments 
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We have developed a new in situ method to calibrate optical tweezers experiments and simultaneously mea- 
sure the size of the trapped particle or the viscosity of the surrounding fluid. The positional fluctuations of the 
trapped particle are recorded with a high-bandwidth photodetector. Next, we compute the mean-square displace- 
ment, as well as the velocity autocorrelation function of the sphere and compare it to the theory of Brownian 
motion including hydrodynamic memory effects. A careful measurement and analysis of the time scales char- 
acterizing the dynamics of the harmonically bound sphere fluctuating in a viscous medium then directly yields 
all relevant parameters. Finally, we test the method for different optical trap strengths, with different bead sizes 
and in different fluids, and we find excellent agreement with the values provided by the manufacturers. The pro- 
posed approach overcomes the most commonly encountered limitations in precision when analyzing the power 
spectrum of position fluctuations in the region around the corner frequency. These low frequencies are usually 
prone to errors due to drift, limitations in the detection and trap linearity as well as short acquisition times re- 
sulting in poor statistics. Furthermore, the strategy can be generalized to Brownian motion in more complex 
environments, provided the adequate theories are available. 



I. INTRODUCTION 

Single Gaussian beam optical tweezers are typically used 
to hold and manipulate a single micron-sized sphere in a 
fluid 0J[3]. The sphere fluctuates by thermal excitations 
within the confines of the harmonic optical trapping poten- 
tial and is employed as a probe of its local environment J4]- 
|7). The positional fluctuations are usually recorded by a pho- 
todiode, which registers the interference pattern between the 
trapping laser scattered by the sphere and the unscattered 
light 1 8 9 ] . Relying on the linearity of the detector response 
ifTUl . the position signal measured in Volts is converted into 
distances by fitting the data to the theory of Brownian mo- 
tion ITiTFFTI . In its simplest form the Einstein-Ornstein- 
Uhlenbeck model describes the thermal fluctuations of a 
sphere in a harmonic potential, where the interaction with 
the surrouding fluid is modelled as an instantaneous Stokes 
friction. In frequency space, the power spectrum of po- 
sitional fluctuations is then represented by an overdamped 
harmonic oscillator with a known viscous drag 7 = Qir-qR, 
where r] is the fluid's viscosity and R the known bead radius. 
Fitting the power spectral density (PSD) yields the Volts-to- 
meter conversion factor, (3, and the trap's spring constant K 
via the corner frequency f c of the Lorentzian. This straight- 
forward and popular calibration method has been extended 
in the past years, to allow also for data calibration with 
spheres of unknown sizes or experiencing an unknown vis- 
cous drag lfT4l - [T9l . The common strategy has been to add 
to the intrinsic thermal fluctuations of the bead, an oscilla- 
tion of known amplitude and frequency by actively moving 
the fluid or the trap relative to the bead, employing a piezo- 
stage or acousto-optic deflectors respectively. 
With the improvement of detection bandwidth |20|, also 
short timescales of optically confined Brownian motion have 
become accessible PTU231 . In principle, three timescales 
characterize the Brownian motion of a bead held in a har- 



monic potential. The longest corresponds to the time dur- 
ing which the particle moves back towards the trap center, 
tk = 6itr)R/K = l/2tvf c . The Brownian particle displaces 
the fluid mass rrif = 47r / Ofi? 3 /3, with pf the fluid density. 
The coupling between the particle and the fluid gives rise 
to the hydrodynamic memory effect created by a fluid vor- 
tex, which builds up around the particle and diffuses over 
the distance of one bead radius within the time Tf = pfi? 2 /r/. 
A third timescale t p = M/6nr)R is obtained by comparing 
the stationary Stokes friction with the inertia of the particles. 
Here, the effective mass M = m p + mf/2 results from a 
calculation including hydrodynamic memory effects. 
We propose to exploit the nature of hydrodynamic inter- 
actions to calibrate high-bandwidth optical tweezers ex- 
periments beyond the Einstein-Ornstein-Uhlenbeck theory 
lfT3l I24H261 . The characteristic timescales of Brownian mo- 
tion, which are exclusively determined by R, r/, K, pf, and 
the particle density p p , can be directly inferred from the ve- 
locity correlation function (VACF) and the mean-square dis- 
placement (MSD) of the bead's fluctuations 1271 . The pre- 
sented method hence does not require a fit to the data rather 
it relies on reading off characteristic features of the measured 
curves. 

In the following, we first give an overview of the experi- 
mental requirements. Then, we demonstrate that the MSD 
and VACF are indeed only dependent on the characteristic 
timescales of confined Brownian motion. Next, we detail 
the calibration procedure on a representative dataset, and fi- 
nally test the method for experiments performed with differ- 
ent trap stiffnesses, different bead sizes and different fluids. 



H. MATERIALS AND METHODS 

The minimal requirements on instrumentation l28l for 
quantitative optical trapping experiments are: (i) a stable 
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light microscope equipped with a high numerical aperture 
(NA) objective for simultaneous optical visualization and 
trapping of the probing sphere, (ii) a highly stable trapping 
laser with an expanded beam for slight overillumination 
of the objective's back aperture, (iii) a photodiode opti- 
mally sensitive to the wavelength of the trapping laser for 
high-resolution and time -resolved position detection 1201 , 
(iv) high-bandwidth signal amplification, and (v) a data 
acquisition board that allows for a sampling rate of at least 
1 MHz. 

It is also convenient to precisely control by, for example a 
three-dimensional (3D) piezo stage, the distance between 
the sphere under study and boundaries in the sample 
chamber lfT0ll29l . 

More specifically, our instrument includes a diode-pumped, 
Gaussian beam, ultra-low noise Nd:YAG laser with a 
wavelength of A=1064 nm, and a maximal light power of 
500 mW in continuous wave mode (IRCL-500-1064-S, 
CrystaLaser, USA). Best trapping efficiency is achieved 
by expanding the effective laser beam diameter 10 times 
with a telecentric lens system (Sill Optics, Germany). The 
IR-beam is reflected by a dichroic mirror (AHF Analysen- 
technik AG, Germany) into the high numerical aperture 
(NA=1.2) of a 60x water-immersion objective (UPLapo/IR, 
Olympus, Japan), which focuses the laser down to its 
diffraction limit into the object plane of the microscope and 
creates the optical trap. The resulting beam waist, Wq, is ap- 
proximately 118 nm in water with a refractive index n=1.32 
at A=1064 nm ll30l . The choice of a water-immersion 
objective lens offers a longer working distance of up to 280 
urn compared to oil-immersion lenses, and guarantees a 
stable uniform trap through the entire sample chamber. The 
sample is mounted onto an xyz-piezo scanning table (P-561, 
Physikalische Instmmente, Germany) for manipulation and 
positioning. The piezo-stage with controller (E-710 Digital 
PZT Controller, Physikalische Instmmente, Germany) has 
a travel range of 100 um in 3D with a precision of ~1 nm. 
The laser light focused by the objective lens is collected 
with a second objective lens (63X, Achroplan, NA = 0.9, 
water-immersion, Zeiss, Germany), and focussed onto an 
InGaAs quadrant photodiode with an active area of 2 mm in 
diameter (QPD, G6849, Hamamatsu Photonics, Japan). The 
QPD is fixed to a x-y translation stage (OWIS, Germany) 
for manual centering of the detector relative to the IR-beam. 
The QPD signals are fed into a pre-amplifier (Offner MSR- 
Technik, Germany), which provides two differential signals 
between the segments and one signal that is proportional 
to the total light intensity. This allows the detection of 
the particle's position in 3D. Pre-amplification of the QPD 
signals at 20 V/mA with 0.67 A/W photosensitivity leads 
to a voltage of 13.4 V/mW. Subsequently, differential 
amplifiers (Qffner MSR-Technik, Germany) adjust the 
pre-amplifier signals for optimal digitalization by the data 
acquisition board with a dynamic range of 12 bits (NI-61 10, 
National Instruments, USA). In accordance with the total 
laser intensity impinging on the QPD, amplification of the 
QPD signal is chosen to span the maximal dynamic range 
of the acquisition card. The amplifier, with a maximal gain 



of 500, has a cut-off frequency around 1 MHz. 
The sample chamber consists of a custom-made flow cell. 
A coverslip (thickness ~ 130 um) is glued to a standard mi- 
croscope slide by two pieces of double-sided tape arranged 
to form a ~ 5 mm wide and ~ 80 um thick channel with a 
volume of ~ 30 pi. After loading with a dilute suspension of 
microspheres, the flow-cell is mounted upside down on the 
3D piezo-stage of our custom-made inverted microscope. 
To achieve the high temporal resolution required for the 
present calibration method, it is crucial to choose the 
detector material to be optimally sensitive to the wavelength 
of the trapping laser 11201 , which in our case is InGaAs. 
Furthermore, a smaller active area results in a lower dark 
current and a higher signal-to-noise ratio at short times. 
For 3D position detection, we chose a QPD diameter of 2 
mm onto which we focussed the trapping laser spot. An 
optimally focussed laser spot yielded a temporal resolution 
close to 1 us. If 3D position detection is not crucial, 
even higher temporal resolution down to 10 ns in a single 
spatial direction can be achieved using balanced amplified 
photodetectors with two well-matched photodiodes with a 
diameter of 0.3 mm each, and an ultra-low noise, high-speed 
transimpedance amplifier [22 1 . 



III. BROWNIAN MOTION 

The individual trajectory of a particle undergoing Brown- 
ian motion is in general unpredictable. Nevertheless, corre- 
lation functions such as the velocity autocorrelation function 
(v(t)v(Q)) and the mean-square displacement (Ax 2 (t)) = 
([x(t) — x(0)] 2 ) are theoretically accessible and describe 
the statistical properties of the dynamics. Both functions 
are connected via (u(t)u(O)) = (l/2)d 2 (Aa; 2 (£))/di 2 . The 
displacement of the particle relative to the trap center as a 
function of time, x(t), obeys the Langevin equation 

m p x{t) = F th {t) + F fl {t) - Kx{t) , (1) 

where m p = 47rp p i? 3 /3 corresponds to the particle's mass, 
Ffa denotes the thermal random noise and Ff r the friction 
forces acting on the particle. The friction force, known as 
Basset force, consists of Stokes' law and additional terms 
accounting for the fluid vortex developing around the parti- 
cle 

2 

Ffr(i) = - 6irqRx{t) - -TrR 3 p f x(t) 

and is valid for spherical particles with no-slip boundary 
conditions on their surface lUTl . Here, the last contribution 
depends on the entire history of the particle's motion and 
therefore it is referred to as "hydrodynamic memory". 
The thermal force describes the random collisions of the 
surrounding solvent molecules with the particle. While its 
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mean vanishes by symmetry arguments, its autocorrelation 
is given by 



(F th (t)F th (t / )) =k B T 7 (\t-t , \), 



(3) 



where j(t) is a kernel describing the retarded friction 
forces. The explicit form of 7(f) requires the use of dis- 
tributions, but its Laplace transform (convention x(s) — 
/ °° x(t) exp(-st)dt) takes the following form |23l l3"fl |32| 

■y(s) = 6nrjR[l + y/sTi] + snif/2 . (4) 



STfJ 

Here, the constant describes the delta-correlated white noise 
usually assumed for Brownian motion. The last term can be 
intrepreted as added mass contibuting to the effective mass 
M, whereas the non-analytic square-root accounts for the 
hydrodynamic memory. 

From the Langevin equation one derives the VACF in the 
Laplace domain, normalized to its initial value (v(Q)v(0)) = 
k B T/M, as 



(v{s)v(0)) _ s 
(v(0)v(0)) ~ W) 



(5) 



with the polynomial P(s) = s 2 + s(l + y / sTf)/T p + l/(r p Tk). 
The inverse Laplace transform is achieved observing that 
the Laplace transform of exp(z 2 t / Yf)erfc(z Wt / Yf) is l/(s+ 
z^Js/Tf). The VACF can be evaluated by partial fraction de- 
composition of Eq. |5]) 



where A« 



(v(s)v(0)) 
(v(O)v(O)) 

■■ lim./:_._ 




(6) 



s(s + Zj ^PFf)/P{s) 



rifc^j [ z j/( z j ~ z k)] ar, d the Zj are the four roots of the poly- 
nomial P(s) — Ylj =1 (\/s + Zj/s/if). Finally, the normal- 
ized VACF is given by 



M*MQ)) 

(v(0)v(0)) 



= Aj exp(z 2 i /rf)erfc(zj 



(7) 



Since in the Laplace domain, (v(s)v(O)) = s 2 (Ai 2 (s))/2, 
a similar partial fraction decomposition may be employed to 
calculate the MSD. After normalization to its plateau value, 
(Ax 2 (t -> 00)) = 2k B T/K, the MSD reads 



(Ax 2 (t 
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^(A i /4)exp(z 2 VTf)erfc(z J -v / V^) • (8) 



Note that Eqs. (j7| and ([8| only depend on the ratio t/rf. The 
trap strength, and the mass of the particle are encoded in the 
dimensionless roots of the polynomial P(s) and in the am- 
plitudes Aj . 

The above description is valid in the case of a spherical parti- 
cle performing Brownian motion in a simple fluid in infinite 



space [26|. Additional time scales arise in more complex 
systems. For example, when approaching a surface, the time 
scale t w = Pfh 2 /i], with h the distance to the surface, de- 
scribes the influence of the wall on the Brownian Motion. 
ET1 l29l . Furthermore, a viscoelastic fluid has a timescale 
that marks the transition from its purely viscous to elastic 
regime (33). However, in such cases, no analytical expres- 
sions for the MSD and VACF are available and their calcu- 
lation has to rely on numerical Fourier transforms. In the 
following section, we show how Brownian motion in an op- 
tical trap may be used to calibrate the optical tweezers only 
in the case where analytic expressions exist. 



IV. PROCEDURE 

Interestingly, as highlighted by Eqs. |7]i and ([8]), the 
shapes of the normalized MSD and VACF are only deter- 
mined by the ratios T p /rf = (2p p /9pf + 1/9) and tr/V? = 
6ttt] 2 / pfKR. Since the temperature and densities of the 
fluid and the bead are usually known or can be easily deter- 
mined by other means, the shape of the curves are sufficient 
to determine the ratio tk/t{ and the timescale Tf, Therefore, 
from reading out r p , Tf, tk, the scale factor (3 can be inferred. 
Eventually, the trap stiffness K and the bead radius R or the 
viscosity r\ follow from the definition of Tf and r K . 
Fig. 1 displays the uncalibrated (Ax 2 (t)) v and (v(t)v(0)) v 
calculated directly from the position signal of an optically 
trapped sphere in water measured in Volts during 20 s at a 
sampling frequency fg = 1 MHz. Assuming a linear detec- 
tor response ifTUl , with ft^ 1 the detector sensitivity in V/m, 
the actual position of the particle is given by 



Xj {t) = px)(t), j = l...N. 



(9) 



where x y is the voltage measured by the detector for N ac- 
quired data points (here N = 2 • 10 7 ). The sampling fre- 
quency determines the time interval between two data points 
St = 1/ fs = 1 |JS. The mean-square displacement is calcu- 
lated from the raw data as: 



1 N 

( Ax2 w) v = ^E[ a;V ( < +^)- a;V ^) 



(10) 



3=1 



and the VACF is determined from the increments of the mea- 
sured voltages through 

1 N 

(v(t)v(0)) v = ]T [x^t + ti + St) x"(t + tj ) 
3=1 

x x y {t 3 +St) -x J {tAl) 

For calibration, two values can be directly read off the 
raw MSD data (Fig. 1(a)): first, the long-time plateau 
(Aa: 2 (oo)) v , which allows the normalization of the data 
(inset), and second, the time tf/ 2 at which (Ax 2 (t)} y 
reaches half of its maximum. The MSD is essentially the 
Fourier transform of the PSD, most commonly used for 
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the calibration of optical traps in the frequency domain 
lfTTHT3l . Similarly, the VACF graph displays two readily 
accessible characteristic features (Fig. 1(b)): the time to at 
which it crosses zero and becomes negative, and the value 
of its minimum (v(t m i n )v(0)) . The anticorrelations in the 
VACF originate from the harmonic restoring force of the 
optical trap and are discussed later in Fig. 3(b). The earlier 
they occur, the stronger the trap. For weaker traps, a third 
interesting feature appears in the positive short-time values 
of the VACF, which decays with a t~ 3 / 2 power-law (Fig. 
1(b), inset). The amplitude of this well-known long-time tail 
(LTT) is given by B = k B J 'Jpf/^-n-qf/ 2 and depends 
only on the fluid properties GT1I251 . 

The ratio t 1 / 2 /to determined from the experimental data 
uniquely sets t^/ts. However, at this step, a comparison 
between the theory and the data in Fig. 1 requires fitting 
with at least 3 variables, f3, K and R or ?/. To reduce 
the number of fit parameters, we compute the theoretical 
functions given by Eqs. (|7]i and ([8]l for a wide range of 
ratios of T^/rf, but fixed r p /rf, hence a given particle and 
fluid density. 

The theoretically possible values for t 1 / 2 as a function of 
Tk/rf are then determined from (Ax 2 (ti/ 2 ))/(Ax 2 (t — > 
oo)) = 1/2, and t is given by (v(t )v(0))/(v(0)v(0)) = 0. 
Consequently, the value of rk/Vf corresponding to an 
experimentally measured value of ti/ 2 /to can be directly 
read off from the bottom graph in Fig. 1(c). Once the ratio 
7\/Tf is known, the values for t\i 2 jri and to/ r f are taken 
from the respective graphs on top in Fig. 1(c), allowing 
the inference of Tf, and subsequently r k , K and either R or 77. 

Finally the conversion factor (3 follows by either: (i) over- 
laying the theoretical and experimental plateau values of the 
MSD: /3 2 SD = 2k B T/(K (Ax 2 (oo)) v ), (ii) matching the 
minimum value (v(t min )v(0)) of the experimental data to 
the minimum of the normalized theoretical curve: /3vacf = 
(v{t min )v(0)) / (v(t min )v(0)) y , or (iii) matching the ampli- 
tude of the LTT measured in V/s 1 / 2 to its theoretical value: 

Applying the described calibration procedure to the exem- 
ple dataset presented in Fig. 1 results in the following steps: 
ty 2 /t = 8.1 ->• TK/Vf = 219 (Fig. 1(c), bottom) ->■ 
t / Tf = 19.7 and t 1/2 /r f = 162 (Fig. 1(c), top) 
T f = 2 us -» r k = 440 us — > for water with r\ = 0.98 cP 
at T = 21° C, we obtain R = yV"f//°f = L4um and 
K = 6n P fR 3 /(T { T k ) = 59uN/m = 0.059 pN/nm, yielding 
three values for the calibration factor f3: /?msd = 9.4nm/V, 
A/acf = 8.8nm/V, and /3 LLT = 8.8nm/V. A fit to the PSD 
as discussed in ref. [ 1 3 1 results in /3pso = 9.8nm/V (data 
not shown), in good agreement with the three values com- 
puted by the method presented here. However, a reliable fit 
of the PSD can only be achieved with two unknown parame- 
ters, typically K and beta. Reading off characteristic features 
from the data rather than fitting to the data allows to deter- 
mine a third parameter. Possible discrepancies between the 
three, differently obtained /3-values and /?psd point to the 



limits in the reading accuracy of each characteristic feature, 
as well as instrumental limitations such as mechanical drift 
and electronic noise. They are discussed in the next sec- 
tion on measurements with different trap stiffnesses, radii 
and viscosities. 



V. TESTS 

Changing the trap stiffness: Fig. 2 displays the calibrated 
MSDs and VACFs of the same resin sphere immersed in wa- 
ter, for various K. The trap stiffness was increased by in- 
creasing the incoming trapping laser power. The black lines 
give the corresponding theoretical predictions. As expected, 
all curves display the same long-time tail, and the plateau 
of the MSD as well as the zero-crossing in the VACF ap- 
pear the earlier given a stronger trap and hence shorter r k . 
Also the minimum of the VACF is deeper, the higher K. In 
contrast, the LTT, which is influenced by the fluid's density 
and viscosity, disappears from our measurement window, as 
motion becomes dominated by the trap. As a consequence, 
depending on the experimental conditions, the features of the 
MSD and VACF may be more or less pronounced. In gen- 
eral, the plateau in the MSD and ti/ 2 can be best determined 
the earlier they are reached. In weaker traps, (Ax 2 (oo)} 
occurs later and may be noisier due to mechanical drift. For 
optimal readability, the zero-crossing in the VACF has to oc- 
cur at to > 2 us, above the resolution limit of our detector. 
Furthermore, (v(t m i n )v(0)) y has to be pronounced enough 
to be distinguishable from the electronic noise floor, which 
here is around 10 6 nm 2 /us 2 . Such conditions are fulfilled for 
1 uN/m < JlT < 500 uN/m, the typical stiffness range of op- 
tical traps. When r k > 500 us a clearly visible long-time tail 
exists and _B V can in turn be evaluated accurately. The pa- 
rameters K and the detector sensitivity obtained by the 
present method are listed in Fig. 2(c). The highest sensitiv- 
ity and the best match between /?msd and /3y^ CF is obtained 
for strong traps, when all respective features are easily dis- 
tinguishable. 

Changing the particle size: Next, we decreased R and 
adapted the laser power in order to have approximately simi- 
lar trap stiffnesses for each bead (Fig. 3). Changing the bead 
size influences the measured signal in an non-linear way and 
at several levels. First, and most straightforwardly, hydro- 
dynamic backflow is enhanced with R 2 via r t , while con- 
finement occurs earlier for bigger particles through r k . Sec- 
ond, the ratio between bead size and laser wavelength, very 
strongly influences the optical forces acting on the particle. 
Empirically and according to theory, forces are strongest for 
2R ~ X/n ~ 800 nm BU . Third, the ratio R/w deter- 
mines the total light intensity impinging on the detector, as 
well as the interference pattern between scattered and non- 
scattered light. These parameters influence the detector sen- 
sitivity /3 _1 , as well as the noise floor of the signal, which 
appears to be higher for smaller beads (Fig. 3). Apart from 
determining K and /3 _1 , we also inferred R directly from 
the data. All obtained parameters are given in the table of 
Fig. 3, showing that the detector sensitivity is reduced for 
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FIG. 1. (Color online) (a) Lin-log representation of the uncalibrated MSD of an optically trapped melamine resin microsphere in water 
(R = 1.4 urn, p p = 1510 g/1, pi = 1000 g/1, r p /r f = 0.45) yielding (Aa; 2 (oo)) v = 1.54 mV 2 and t 1/2 = 327 jis. Inset: Log-log 
representation of the same MSD normalized by its plateau value, (b) Lin-log representation of the corresponding uncalibrated VACF, with 
t = 39.3 us and (v (t mln )v(0)) v = -2.64 ■ 10 6 mV 2 /s 2 . Inset: Log-log representation of \(v(t)v(0)} v \. The amplitude of the t' 3/2 
long-time tail reads B v = 0.65mV 2 /s 1 / 2 . Data were blocked in 10 bins per decade. Error bars give the standard error on the mean from 
blocking, (c) Top: Theoretical ratios ti/2/rf and to/n as a function of Tt/rt for a resin sphere (p p = 1510 g/1) in water {pi — 1000 g/1). 
The trap relaxation time increases from left to right, while the trap stiffness decreases. Bottom: Corresponding theoretical values of ti/2/to. 
Each value necessary to calibrate the data shown as an example is highlighted by an arrow. 
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FIG. 2. (Color online) Calibrated MSD (a) and VACF (b) of single resin sphere (p p = 1510 g/1, R = 1.4 urn) in water (p f = 1000 g/1, 
r\ = 0.98 cP, T = 21° C, Tp/rt = 0.45, Tt = 2 us) for increasing K. Data were blocked in 10 bins per decade. Error bars give the standard 
error on the mean from blocking, (c) Overview of the parameters obtained from calibrating the data. 



2R < A. The obtained i?g t -values for an individual parti- 
cle are in excellent agreement with the bead sizes, R pro- 
vided by the manufacturer. However, for bead radii yielding 
Tf = pfR 2 /rj < 1 us, hydrodynamic effects become less 
prominent in our measurement window, which impairs the 
reading off accuracy of the short-time features, £? v and to. 
A simultaneous quantification of R and K then becomes a 
challenge. Here, a detector with an improved temporal res- 
olution as used in E2l would allow to measure even smaller 
particle sizes. 

Changing the fluid viscosity: Finally, we tested our 
method to calibrate the optical trap and simultaneously mea- 
sure the viscosity of a fluid, when the bead size is known 
and the laser power is held constant. Different fluids also 



alter the measured signal in a non-linear way. Firstly, the re- 
fractive index and the absorption coefficient of the medium 
influence the trapping forces and the light pattern reaching 
the detector. Secondly, a higher viscosity will yield weaker 
trap stiffnesses as quantified by r^. Finally, the shape of the 
VACF (Eq. 0) and the MSD (Eq.Q) are most dependant 
on rj 2 through T^/Vf. In contrast, a shift of r p /rf correspond- 
ing to a change in the fluid density has a negligible effect on 
the MSD and VACF. Fig. 4 displays the calibrated MSD (a) 
and VACF (b) of a resin sphere (R — 1.47 um) immersed ei- 
ther (i) in the low viscosity solvent acetone, (ii) in water, or 
(iii) in 30% glycerol, which is 2 to 3 times more viscous than 
water. The trap stiffness K, the viscosity r\ and the conver- 
sion factor were determined from the presented method 
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FIG. 3. (Color online) Calibrated MSD (a) and VACF (b) of single resin particles (p p — 1510g/l) with different sizes in water (pf = 1000g/l, 
r\ = 0.98 cP, Tp/Vf = 0.45). Data were blocked in 10 bins per decade. Error bars give the standard error on the mean from blocking, (c) 
Overview of the parameters obtained from calibrating the data. 
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FIG. 4. (Color online) Calibrated MSD (a) and VACF (b) of resin particles in acetone (pf = 790 g/1), water (pf = 1000 g/1) and glycerol 
(pf — 1066 g/1). Data were blocked in 10 bins per decade. Error bars give the standard error on the mean from blocking, (c) Overview of 
the parameters obtained from calibrating the data. 



and are listed in Fig. 4(c). The obtained fluid viscosities are 
in good agreement with viscosimetry measurements, which 
valued rj — 0.32 cP for aceton, rj — 0.98 cP for water, and 
rj = 2.3 cP for 30% glycerol in water at ambient tempera- 
ture. 



VI. DISCUSSION AND CONCLUSION 

In summary, we have established a procedure to calibrate 
optical trapping experiments with a micron-sized particle 
immersed in a viscous fluid. The method allows the determi- 
nation of the position detector sensitivity, the trap stiffness 
and simultaneously the particle size or the fluid viscosity 
with high precision. When neither 77 nor R are known, we 
suggest performing two consecutive measurements at two 
different trap strengths. For a weak trap, rj can be extracted 
from the short-time hydrodynamic features given by the am- 



plitude B of the t~ 3 / 2 power-law. Upon knowing the fluid 
properties, the procedure can be pursued as described. Fur- 
thermore, the presented method can be applied equivalently 
in both directions, x and y, perpendicular to the optical axis, 
yielding K x , K y , f3 x and f3 y . For the dimension, z along the 
optical axis, the signal-to-noise ratio, and hence bandwidth 
of the position detector is usually not sufficient to measure a 
smooth VACF exhibiting a clearly detectable zero-crossing 
(see supplemenary information in |23|). Nevertheless, af- 
ter having determined the bead radius and/or fluid viscosity 
through analysis of the x or y dimension, (3 Z and K z can 
be straightforwardly inferred from the plateau value of the 
MSD and the corner frequency f c = 1/2ttt^ of the PSD 
lfT3l . which is much less prone to high-frequency noise than 
the VACF. Measuring in 3D should allow the assesment of 
anisotropics in the bead sizes as well as in the surrounding 
medium. The proposed approach should be applicable in mi- 
crorheology for the study of viscoelastic fluids, in which at 



7 



least one additional timescale marks the transition between 
purely viscous and elastic behavior 11331 . Also, when the 
Brownian particle approaches a surface, a new characteristic 
time can be extracted from the data, indicating the distance 
between the particle and the surface ETI . 
For example, when the Brownian particle approaches a sur- 
face, the time needed by the hydrodynamic vortex to diffuse 
from the particle to the surface can be read from the data. 
Measuring this time yields the distance between the particle 
and the surface between a purely viscous and an elastic be- 



havior of the fluid ll33ll . and characterizes Brownian motion. 
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